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Abstract 

We conjecture that the asymptotic behavior of the numbers of solid 
(three-dimensional) partitions is identical to the asymptotics of the three- 
dimensional MacMahon numbers. Evidence is provided by an exact enu- 
meration of solid partitions of all integers < 68 whose numbers are re- 
produced with surprising accuracy using the asymptotic formula (with 
one free parameter) and better accuracy on increasing the number of free 
parameters. We also conjecture that similar behavior holds for higher- 
dimensional partitions and provide some preliminary evidence for four and 
five-dimensional partitions. 
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The purpose of computation is insight, not numbers. - Richard Hamming 

1 Introduction 

Partitions of integers appear in large number of areas such number theory, combi- 
natorics, statistical physics and string theory. Several properties of partitions, in 
particular, its asymptotics (the Hardy- Ramanuj an- Rademacher formula) can be 
derived due to its connection with the Dedekind eta function which is a modular 
form [HE]- In 1916, MacMahon introduced higher-dimensional partitions as a 
natural generalization of the usual partitions of integers [3]. He also conjectured 
generating functions for these partitions and was able to prove that his generating 
function for plane (two-dimensional) partitions was the correct one. However it 
turned out that his generating function for dimensions greater than two turned 
out to be incorrect. Even for plane partitions, one no longer has nice modular 
properties for the generating function. Nevertheless, the existence of a generating 
function enables one to derive asymptotic formulae for the numbers of plane par- 
titions The inability to do the same with higher-dimensional partitions (for 
dimensions > 2) has meant that these objects have not been studied extensively. 
The last detailed study, to the best of our knowledge, is due to Atkin et. al. [5]. 

Higher-dimensional partitions do appear in several areas of physics (as well 
as mathematics) and thus it is indeed of interest to understand them better. It 
is known that the infinite state Potts model in [d + 1) dimensions gets related 
to d- dimensional partitions P,[7j. They also appear in the study of directed 
compact lattice animals [H]; in the counting of BPS states in string theory and 
supersymmetric field theory ^^U\ . For instance, it is known that the numbers of 
mesonic and baryonic gauge invariant operators in some = 1 supersymmetric 
field theories get mapped to higher- dimensional partitions [9]. The Gopakumar- 
Vafa (Donaldson-Thomas) invariants (in particular, the zero-brane contributions) 
are also related to deformed versions of higher- dimensional partitions (usually 
plane partitions) [TTl[T2](see also [T3]). 

In this paper, we address the issue of asymptotics of higher-dimensional parti- 
tions as well as explicit enumeration of higher-dimensional partitions. Our work 
builds on the seminal work of Mustonen and Rajesh on the asymptotics of solid 
partitions [H]. lack of a simple formula for the generating functions of these 
partitions has been a significant hurdle in their study. The conjectures on the 
asymptotics of higher dimensional partitions given in this paper, even if partly 
true, would constitute progress in the study of higher- dimensional partitions. The 
conjecture on the asymptotics was arrived upon serendipitously by us when we 
found that a one-parameter formula for solid-partitions derived using MacMa- 
hon's generating function worked a lot better than it should. To be precise, a 
formula that was meant to obtain an order of magnitude estimate (for solid par- 
titions of integers in the range [50,62]) was not only getting the right order of 
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magnitude but was also correct to 0.1 — 0.5% (around 3-4 digits). The main 
conjecture discussed in section 3 is a natural outgrowth of this observation. The 
exact enumeration of solid partitions was possible due to an observation that lead 
to a gain of the order of 10"^ to 10^ enabling us to exactly generate numbers of 
the order of 10^^ — 10^^ in reasonable time. 

The paper is organized as follows. Following the introductory section, section 
2 provides the background to problem of interest as well as fixes the notation. 
Section 3 deals with asymptotics of higher- dimensional partitions. This done by 
means of two conjectures. We provide some evidence towards these conjectures 
with a fairly detailed study of solid partitions using a combination of exact enu- 
meration as well as fits to the data. Section 4 provides the theoretical background 
to the method used for the exact enumeration of higher-dimensional partitions. 
We conclude in section 5 with some remarks on extensions of this work. In ap- 
pendices A we work out the asymptotics of MacMahon numbers. Appendix B 
provides an 'exact' asymptotic formula for three-dimensional MacMahon num- 
bers. In appendix C we present several tables that includes our results from exact 
enumeration as well some details of the fits for solid partitions. 

2 Background 

A partition of an integer n, is a weakly decreasing sequence (oq, Oi, 02, . . .) such 
that 

• Z)i «i = ^nd 

• ai+i < ai V i. 

For instance, (2, 1, 1) is a partition of 4. Define pi{n) to be the number of parti- 
tions of n. For instance, 

4 = 3 + 1 = 2 + 2 = 2 + 1 + 1 = 1 + 1 + 1 + 1 =^ pi(4) = 5. (2.1) 

A slightly more formal way definition of a partition is as a map from Z>o to 
Z>o satisfying the two conditions mentioned above. This definition enables one 

to generalise to higher dimensional partitions. A d- dimensional partition of n is 
defined to be a map from Z>q to Z>o such that it is weakly decreasing along all 
directions and the sum of all its entries add to n. Let us denote the partition by 
{o,h,i2,...,id)- The weakly decreasing condition along the r-th direction implies that 

0'ii,i2,...,ir+l,-,id - <^ii,i2,-,ir,-id ^ {h,i2,---, id) ■ (2-2) 

Two-dimensional partitions are also called plane partitions while three-dimensional 
partitions are also called solid partitions. Plane partitions can thus be written 
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out as a two-dimensional array of numbers, aij. For instance, the two-dimensional 
partitions of 4 are 



2 111 11 11 



3 1 f 2 2 2 2 1 1 2 1 I mi 111 11 1 (2.3) 



1 ^ 1 , 



Thus we see that there are 13 two-dimensional partitions of 4. Let us denote 
by pd{n) the number of d- dimensional partitions of It is useful to define the 
generating function of these partitions by (pd(0) = 1) 

oo 

P,(g)^^Mn)g". (2.4) 

n=0 

The generating functions of one and two-dimensional partitions have very nice 
product representations. One has the Euler formula for the generating function 
of partitions 

Piiq) = ffso^. ^ , (2.5) 

nn=i(i - ? ) 

and the MacMahon formula for the generating function of plane partitions 

^2(5) = fr^T^ — ^ • (2.6) 

nn=i(i - 1) 

MacMahon also guessed a product formula for the generating functions for d > 2 
that turned out to be wrong [5j. His guess is of the form 

^ 00 

We will refer to the numbers md{n) as the (i-dimensional MacMahon numbers. It 
is easy to see that Mi(g) = Pi{q) and M2(g) = P2{<l)- However Md{q) 7^ Pd^q) for 
d > 2. An explicit formula (given by Atkin et. al. [5J or the book by Andrews |15] ) 
for the number of (i-dimensional partitions of 6 is 

p.(6) = l + 10d+27(^]+28Q+llQ + (^] . (2.8) 



Then, one can show that 



mdi6)-Pd{6)=(f] + (f] , (2.9) 



3/ V4 

which is non- vanishing for d > 3. Thus the MacMahon generating function fails 
to generate numbers of partitions when d > 3. 



^We caution the reader that there is another definition of dimensionaUty of a partition that 
differs from ours. For instance, plane partitions would be three-dimensional partitions in the 
nomenclature used in Atkin et. al. [5| while we refer to them as two-dimensional partitions. 
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2.1 Presentations of higher-dimensional partitions 

There are several ways to depict higher dimensional partitions. Recall that there 
is a one to one correspondence between (one-dimensional) partitions of n and Fer- 
rers (or Young) diagrams. The partition of 4 corresponding to 3 + 1 corresponds 
to the Ferrers diagram 



Similarly, the plane partition \ can be represented by a Young tableau (i.e., 
a Ferrers diagram with numbers in the boxes) or as a 'pile of cubes' stacked 
in three dimensions (one of the corners of the cubes being located at (0,0,0), 
(0, 0, 1), (0, 0, 2) and (1, 0, 0) in a suitably chosen coordinate system) 
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Similarly, (i-dimensional partitions can be represented as a pile of hypercubes in 
{d+1) dimensions. 

We refer the reader to the work by Stanley (and references therein) for an 
introduction to plane partitions [HKT^. The book by Andrews [15] provides a 
nice introduction to higher-dimensional partitions. Further the lectures by Wilf 
on integer partitions [TB] and the notes by Finch on partitions [19J are also good 
starting points to existing literature on the subject. 



3 Asymptotics of higher-dimensional partitions 

In this section, we will discuss the asymptotics of higher- dimensional partitions. 
The absence of an explicit formula for the generating function for d > 2 implies 
that there is no simple way to obtain the asymptotics of such partitions. In this 
regard, an important result due to Bhatia et. al. states that [8] 

lim n"'^/'^"'"^ logprf(n) = (i-dependent constant. (3.1) 

n— )-oo 

Conjecture 3.1 The constant in the above formula is identical to the one for 
the corresponding MacMahon numbers. 

lim n~d^ \ogpd{n) = lim n~'^ \ogmd{n) = — ^ 

(3.2) 

For three-dimensional partitions, this becomes a conjecture of Mustonen and 
Rajesh. Mustonen and Rajesh used Monte-Carlo simulations to compute the 
constant and showed that it is 1.79 ± 0.01 [H]. This is compatible with the 
conjecture since /3J ~ 1.78982. 



d C{d+1) 



d+l 



id) 
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It is important to know the sub-leading behavior of the asymptotics of higher- 
dimensional partitions in order to have quantitative estimate of errors. This is 
something we will provide in the next subsection. Before discussing the asymp- 
totic behavior of the higher- dimensional partitions, it is useful to know the asymp- 
totic behavior of the MacMahon numbers. A calculation shown in appendix A 
gives their sub-leading behavior. One obtains 

d 

log main) ~ ^ Pi"^^ + 7^'^^ log n + S'-'^^ . (3.3) 

r=l 

The constants (Si*^^ and 7*-'^-' have been computed for c? = 3, 4, 5 in appendix A. 



3.1 Towards a stronger conjecture 

The number of d- dimensional partitions of n can be obtained from the generating 
function Pd{q) by inverting Eq. (12. 4p 

Pdin) = ^ £ P.ie^y) e-'^^y dy . (3.4) 

Suppose we knew all the singularities of the function Pd{q)- The integral can 
be then be evaluated (at large n), for instance, by the saddle point method 
and adding up the contribution of all singularities thus obtaining an asymptotic 
formula for pd{n). The singularities are usually obtained by looking at product 
formulae of the form 

00 

^d(g) = n(i-^")""'"^"^ ■ (3-5) 

71=1 

The exponents a^'^\n) can be determined for those values of n for which pd{n) has 
been determined. If all the a^^\n) are positive, then it is easy to see that -Pd(g) is 
singular at all roots of unity - this leads naturally to the circle method of Hardy 
and Ramanujan [T]. However, for d > 2, this turns out to be false. For instance, 
a^^\l5) = —186 is the first exponent that becomes negative for c? = 3 pUi see 
Table 1]. We will assume that the singularities of Pd{(l) continues to occur at 
roots of unity. In particular, we will see that the Bhatia et. al. result implies 
that for large enough n, one has 

aW(n) = 0{n''-^) , (3.6) 

with a'-'^\n) > 0. Let us assume that the dominant term in a saddle point 
computation of the integral in Eq. (13. 4p occurs near q = 1. 

Proposition 3.2 The Laurent expansion 0/ logPd(e~*) in the neighbourhood of 
t = is of the form 

— logPd(e~*) = -I- , , ^'^7^ , , -I- h — -I- non-sinqular as t , (3.7) 

^ ^ df^ {d-l)t'^~^ t ' V y 
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where Ci, . . . ,Cd o'^e some constants. 

Remark: This is precisely the form of the Laurent expansion for logMd(e~*) 
near t = (see Appendix A). 

A saddle point computation of the integral (13 ■4p is carried out by extremizing 
the function 

hg Pd{e-') +nt . 

The extremum, t^,, which is close to t = for large n, obtained using Proposition 
13.21 is given by 

t*=[^j +■■■ (3-8) 

Plugging in the saddle point value, we see that 
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(^dd^ + sub-leading terms . (3.10) 

We thus recover the bound obtained by Bhatia et. al. [8]. Thus, we see that the 
Bhatia et. al. result combined with the assumption that Pd(e~*) is a meromorphic 
function in the neighborhood of t = with a pole of order d implies Proposition 

m 

A more precise saddle point computation enables us to determine sub-leading 
terms as well and we obtain 

d 

\ogpd{n) ~ P^f^ + ^^"'^ log^ + 5^'^^ + ■■■ , (3.11) 

r=l 



where the constants I3l'^\ ^^'^^ and 5*-'^^ are determined by the constants Cr that 
appear in Proposition 13.21 

Conjecture [XT] implies that Cd = d({d+ 1) - this is the leading coefficient in 
the Laurent expansion of logMrf(e~*) near t = 0. This is equivalent to 

where the ellipsis indicates sub-leading terms in the large n limit. We now propose 
a stronger form of conjecture 13. 1[ 

Conjecture 3.3 The asymptotics of the d-dimensional partitions are identical 
to the asymptotics of the MacMahon numbers. 

d 

\ogpd{n) ~ Pf> + 7^'^) logn + ■ • • , (3.13) 

T- = l 



where pl^^ and 7'-'^^ are as in Eq. (13.31) . 
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61 


17474590403967699 


17472947006257293 


62 


used to fit constant 


28518691093388854 


63 


46453074905306481 


46458506464748807 


64 


75522726337662733 


75542021868032878 


65 


122556018966297693 


122606799866017598 


66 


198518226269824763 


198635761249922839 


67 


320988410810838956 


321241075686259326 


68 


518102330350099210 


518619444932991189 



Table 1: Estimates using the asymptotic formula q^{N). The constant in the 
asymptotic formula is fixed by requiring it to give the exact answer for iV = 62 
- the largest known number of solid partitions at the time of the fit. 



It is easy to see that one can have conjectures that are stronger than conjecture 
13.11 but weaker than conjecture 13.31 by requiring fewer coefficients to match with 
Eq. (13. 3p . Conjecture 13.31 implies that the coefficients, Cr {r = l,...,d) in 
the Laurent expansion in Proposition 13.21 are identical to those of logMd{e~*). 
Equivalently, 

P,(e-*) - Mrf(e-*) = 0{l) , (3.14) 

near t = 0. It also implies that at large n, a^'^\n) behaves exactly like the expo- 
nent that appears in the product formula for d-dimensional MacMahon numbers 
in Eq. (EZD, i.e., 



d 

where the ellipsis indicates terms that vanish as n — )■ oo 



3.2 Evidence for the conjecture 

We will provide evidence by explicitly enumerating numbers for the higher- 
dimensional partitions. In particular, we compute all solid partitions for < 68 
and use the formula provided by Eq. (13.131) as a one-parameter function to fit 
known numbers. The advantage of this procedure is that one doesn't need to go 
to enormously large values of n. In Figures 1, 2 and 3, we compare this formula 
implied by conjecture 13.31 for o? = 3, 4, 5 respectively. Since the values of n that 
we consider are not too large, these fits provide weak evidence that three of the 
conjectured numbers i.e., fyf \ fi'^'^ and 7*^*^^ are probably correct. 
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Figure 1: Plot of n~^^^\ogp3{n) for n E [5,68] (red dots). The blue curve is the 
asymptotic formula normalized to give the correct answer for n = 62 and the 
horizontal line is the conjectured value for n — )■ oo. 

3.3 Solid partitions: a detailed study 

The asymptotic expansion of the logarithm of three-dimensional MacMahon num- 
bers is (with ^ = n 



C(-3)^ 



4 



logma^nj g[jc,^4jj ^ + 2[3C(4)]i/2 ^ 8[3C(4)]5/4 ^ 96 ^^ 

(3.16) 

Using the above formula as a guide, we fit the solid partitions to the following 
three formulae involving up to three parameters (a, c): (^ := n + 6) 

gsW g[dU4JJ n + 2[3^(4)]i/2 ^ 8[3C(4)]5/4 ^ gg^^S^i + a 



gi-^^ ^ + 2[3C(4)]V2 ^ 8[3C(4)]5/4 ^ 96 ^ 



Note that the number of free parameters increases from 1 for the function gs to 
2 for rg and to 3 for S3. We obtain a = -1.544, (a, 6) = (-1.530,-0.028) and 
(a, b, c) = (—3.211, 1.689, 0.257) from the three fits. We use the same functions to 
estimate the values of three-dimensional MacMahon numbers for the same range 
of values using a similar fit. We see that the function S3 has worked almost as 
well as it did for the corresponding MacMahon numbers. In particular, the fit 
gives c = 0.25713 which is different from the one given by MacMahon numbers for 
^ -0.041413. For the MacMahon numbers, the fitted value of c = -0.057621 
which is close to the actual number. This suggests that the coefficient of n^/^ 
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n I logp4(n) 



1.6 - 



1.5 - 



1.4 - 




10 20 30 



Figure 2: Plot of n~^^^\ogpi{n) for n E [5,35] (red dots). The blue curve is the 
asymptotic formula normalized to give the correct answer for n = 30 and the 
horizontal line is the conjectured value for n oo. 

may be different from the one given by the MacMahon numbers. For the values 
of n that we have considered, the dominant contributions are due to the first two 
terms as well as the log term. Hence, we consider this as possible evidence for 
Pr = /3r for r = 1,2. For completeness, we provide the numbers obtained by 
carrying out a five-parameter fit using the numbers in the range [60, 68]. The fit 
gives: 

logp3(n) ~ 1.73 n^/^ + 0.83 n^^^ - 0.90 n^^^ - 1.00 logn - 0.22 . (3.17) 

We also observe that if we used a larger range of numbers, say, n G [50, 68] , we 
obtain large numbers (of order ten or greater) for some of the coefficients. This 
reflects the lack of data for large number more than anything else. 

In an attempt at understanding the accuracy of our numbers better, we car- 
ried out a systematic study of an exact asymptotic formula (in the sense of 
Hardy-Ramanujan-Rademacher for partitions) for three-dimensional MacMahon 
numbers using a method due to Almkvist [2T1[22]. These are discussed in Ap- 
pendix B. One writes 

oo 

k=l 

where 4>k{n) are the contributions from various saddle-points with k = 1 being 
the dominant one. For n = 60, we see that 0i(6O) gets the first nine digits right 
while the sum of the first two terms get eleven digits right. We further broke up 
the contribution of (f)i{n) into several terms. The term that we write as (pf^n) 
is the contribution from the singular part of logM3(e~t) at the dominant saddle 
point located near t = 0. We see that 0f''(6O) gets the first five digits right - 
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log p,(n) 




Figure 3: Plot of n'^^^logp^^n) for n E [5,30] (red dots). The blue curve is the 
asymptotic formula normalized to give the correct answer for n = 25 and the 
horizontal line is the conjectured value for n oo. 

somewhat closer to what we have obtained in our estimates for the numbers of 
solid partitions. 

3.4 An unbiased estimate for the leading coeffficient 

In order to provide an unbiased estimate for the leading coefficient of the asymp- 
totic formula using the exact numbers of solid partition^, we use the method of 
Neville tables (albeit with a slight and obvious modification) [23]. Let 

e° = n~^^^ log p^{n) 

~^/3(Vi-'-)/^ + 7(=^)n-=^/Mogn + 5(3V3/4 ^ (g^^g^ 

x=l 

where we have written the asymptotic formula in the second line using the pa- 
rameters defined in Eq. (13. lip . Further, for r > 1, recursively define 



- (n-rV/^ e"-^- 



n-l 



-in- r)i/4 



(3.19) 



Using the conjectured asymptotic formula for P3(n), we can derive asymptotic 
formulae for e^. The have been constructed so that 

1. lim„_!.oo ej^ tends to a constant that equals for all r. The first sub- 
leading term is proportional to n~'^''^^^/^. Thus a plot of vs n~^'''^^'^^^ 
should be a straight line in the asymptotic limit. 



^We thank the anonymous referee for suggesting that we provide an unbiased estimate of 
the leading coefficient and for asking us to look at the methods discussed in ref. [23] , 
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2. As we increase r, the number of parameters that appear in the asymptotic 
formula for decrease. For instance, one sees that /^g^"* drops out for r = 1: 

4 ~ /3f ^ - /3f _ 2^(3)^-3/4 ^ ^ (4^(3) _ 25(3))^-3/4 ^ (3,20) 
and /^g^'', /33^'' drop out for r = 2: 



e: 



2 . , «{3) 



/3(3) +^(3)^-3/4 iQg^^ (_g^(3) ^^(3))^-3/4 ^ (3,2!) 



An estimate for /3{ has been obtained by carrying out two and three-parameter 
fits to the asymptotic formula given in Eq. (I3.20p . We obtain 



en 



1.793 + 2.099^^1/2 two-parameter fit , - 

1.781 + 0.83n-i/2^ 0.924 log n three-parameter fit ^ ' ' 



A four-parameter fit leads to coefficients that are not of order one. We discard 

this fit as we make the natural assumption that all coefficients are of order one or 

" (3) 

smaller. Using the two different fits, we can estimate /3{ is around 1.78 — 1.79. 

^ (3) 

The wide variation that we observe in suggests that we cannot estimate it 
with the available exact numbers. In figure HJ we have plotted e\ vs n"^/^ along 
with the three-parameter fit. We also observe that (see figure |5]) is oscillat- 
ing between [1.77, 1.81] and hence we cannot estimate any further parameters 
using the data. For completness, we have carried out a similar analysis for the 
MacMahon numbers, m^{n) in the range n G [20, 68] and obtain /3f ^ in the range 
[1.77 — 1.78]. We also observe that does not oscillate as it does for solid 
partitions. 

We conclude that an unbiased estimate for /J^^-* is consistent with conjecture 
13.11 However, given the relatively small values of n that we have used, this only 
constitutes weak evidence at best. There is another result due to Widom et. al. 
who studied the asymptotics of (restricted) solid partitions with Ferrers diagrams 
that fit in a four- dimensional box of size 10^ x p [24J as a function of p. They 
observe that the entropy in the thermodynamic limit deviate^ from a formula 
derived from a MacMahon formula for restricted solid partitions. Should we ex- 
pect a similar behavior for unrestricted solid partitions? The deviation observed 
by Widom et. al. is small. If a similar behavior occurs for unrestricted partitions, 
then conjecture 13.11 would be false. We believe that the exact numbers that we 
have used are not large enough to definitively test conjecture 13.11 However, in 
any case, it is important to note that that the functional form of the asympotics 
continues to hold. 



■^The entropy for fixed boundary conditions was found to be 0.145 instead of the conjectured 
value of 0.139. See Eq. (14) in ref. [24]. 
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Figure 4: A plot of 



vs n 



for n G [20, 68] along with a three-parameter fit. 



4 Explicit Enumeration 

In this section, we discuss the explicit enumeration of higher dimensional par- 
titions. The first program to explicitly enumerate higher-dimensional partitions 
is due to Bratley and McKay [22]. However, we do not use their algorithm but 
another one due to Knuth [20]. We start with a few mathematical preliminaries 
in order to understand the Knuth algorithm as well as our parallelization of the 
algorithm. 

4.1 Almost Topological Sequences 

Let P be a set with a partial ordering (given by a relation denoted by -<) and a 
well-ordering (given by a relation denoted by <). Further, let the partial ordering 
be embedded in the well-ordering i..e, x -< y implies x < y. 

Definition 4.1 A sequence X = {xi,X2, ■ ■ ■ ,Xm) containing elements of P is 
called a topological sequence if fWf 

1. For 1 < j < m and x ^ P , x -< Xj implies x = Xi for some i < j ; 

2. If m > 0, there exists x E P such that x < Xm and x ^ Xi, for 1 <i <m. 

Let us call a j-th position in a topological sequence, X, interesting if Xj > Xj+i. 
By definition, the last position of a sequence is considered interesting. The index 
of a topological sequence is defined to be the sum of all j for all interesting 
positions i.e.. 




(4.1) 
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1.82 



1.81 



4.80 



1.79 



1.78 



0.050 0.055 0.060 0.065 0.070 0.075 
Figure 5: A plot of vs n"^/'' for n G [20, 68]. 

Definition 4.2 y4n almost topological sequence is a sequence that satisfies con- 
dition 1 hut not necessarily condition 2. 

Tlius all topological sequences are also almost topological sequences. This defini- 
tion is motivated by the observation that almost topological sequences do occur 
as sub-sequences of topological sequences. 

4.1.1 An example due to Knuth 

Let P denote the set of three-dimensional lattice points i.e., 

P= {(^,j,A;) I 2,j,A; = 0,l,2,3,...} (4.2) 

with the partial ordering (i, j, k) ^ (i', j', k') ii i < i' , and j < f and k < k' . Let 
us choose the well-ordering to be given by the lexicographic ordering i.e., 

(^,j,A;)<(z',/,fc') (4.3) 

if and only if 

i < i' oi {i = i' and j < j) or {i = i' , j = j' and k < k') . 

The depth of a topological sequence is the number of elements in the sequence. 
Consider the topological sequence (of depth 6) 

X = {(0, 0, 0), (0, 0, 1), (0, 0, 2), (1, 0, 0), (0, 1, 0), (0, 0, 3)} 

where we have indicated the interesting positions in boldface. This sequence has 
index 15 = 4 + 5 + 6. 
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4.2 Topological sequences and solid partitions 



Let dm{n) denote the number of topological sequences of the set P = N™" with 
index n. Further, define dm{0) = 1. As before, let Pm{n) denote the number 
of m-dimensional partitions of n. A theorem of Knuth relates these two sets of 
numbers as follows: 

Theorem 4.3 (Knuth [20]) 

n 

Pm{n) = ^dmik) pi{n - k) . (4.4) 

k=0 

Equivalently, the generating function of m-dimensional partitions decomposes 
into a product of the generating function of the numbers of topological sequences 
and the generating function of one-dimensional partitions. 

PM = D^q) P,{q) , (4.5) 

where 

oo 

Dm{q) ■=^dm{n) q"" . 

n=0 

Since topological sequences are much easier to enumerate, Knuth went ahead 
and wrote a program to generate all topological sequences of index < (for 
some fixed A^). This is the program that was the starting point of our exact 
enumeration. 

We list below the topological sequences of index 2 and 3 when P = (we 
have dropped the comma between numbers to reduce the length of the expression) 



Index 2: {(000) (010)} and {(000)(100)} =^ d^{2) 



Index 3: {(000)(001)(010)} ; {(000)(001)(100)} ; {(000)(010)(020)} 



{(000)(010)(100)} ; {(000)(100)(200)} =^ d^iS) 



Thus, we see that Dslq) = 1 + 2g^ + 5g^ + ■ ■ ■ . We also have Pi{q) = 1 + g + 
2g^ + 3g^ + ■ ■ ■ . Thus, we obtain 

P3(g) = Dsiq) giiq) = l + q + 4q' + W + ■■■ 



4.3 Equivalence classes of almost topological sequences 

We say that two sequences X = (xi, X2, ■ ■ ■ , Xm) ~ Y = {yi, 1/2, ... , ym) are re- 
lated if the elements of Y are a permutation of the elements of X. Of course, 
not all permutations of an almost topological sequence lead to another almost 
topological sequence as some of them violate condition 1 in the definition of a 
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topological sequence. However, even after imposing the restriction to permuta- 
tions that lead to other topological sequences, the relation remains an equivalence 
relation. As an example consider the following three sequences in N^: 

{(0,0,0), (0,0,1), (0,0, 2), (1,0,0)} , 

{(0,0,0), (0,0,1), (1,0,0), (0,0, 2)} , (4.6) 
{(0,0,0), (1,0,0), (0,0,1), (0,0, 2)} . 

It is easy to see that these three sequences form a single equivalence class. How- 
ever, the last two are not topological sequences as they violate condition 2 in the 
definition of a topological sequence, and hence are almost topological sequences. 
We thus choose to work with equivalence classes of almost topological sequences. 

Proposition 4.4 The equivalence classes of almost topological sequences of N'^ 
of depth k is in one to one correspondence with {d — 1) -dimensional partitions 
of k. We shall refer to the {d — 1) -dimensional partition as the shape of the 
equivalence class. 

The ((i— l)-dimensional partition is obtained by placing d- dimensional hypercubes 
(of size one) at the points appearing the almost topological sequence. This is 
nothing but the 'piles of cubes' representation of a (d — l)-dimensional partition. 
In this representation, the precise ordering of the points in the almost topological 
sequence is lost and one obtains the same {d — l)-dimensional partition for any 
element in the same equivalence class. Given a (rf — l)-dimensional partition, 
the coordinates of the hypercubes in the 'piles of cubes' representation give the 
elements of the almost topological sequence. For instance, the equivalence class 
in Eq. (14. 6 P has as its shape the following two-dimensional partition of 4: 



3 

T 




When P = N^, the almost topological sequences of P are standard Young 
tableaux. Given an almost topological sequence of with shape A with n boxes, 
the standard Young tableau is obtained by entering the position of the box in 
the almost topological sequenc^. It is easy to see that this map is a bijection. 
It is an interesting and open problem to enumerate the number of almost topo- 
logical sequences given a shape for higher- dimensions. We did this by generating 
all topological partitions of a given index and sorting them out by shape. How- 
ever, this is an overkill if one is interested in enumerating topological sequences 
associated with a particular shape. 

^Recall that a Young tableau is a Ferrers diagram with boxes fihed in with numbers. A 
standard Young tableau has numbers from (1, . . . , n) such that the numbers in the boxes increase 
as one moves down a column or to the right. 
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Depth 12 


14 


15 


17 


Nodes 28680717 


1567344549 


12345147705 


856212871761 


Shapes 1479 


4167 


6879 


18334 



Table 2: Number of equivalence classes at various depths (equal to the number 
of plane partitions) for counting topological partitions of N'^. 

4.4 Programming Aspects 

The explicit enumeration of topological sequences to generate partitions was first 
carried out Knuth who enumerated solid partitions of integers < 28 [20] . This was 
extended to all integers < 50 by Mustonen and Rajesh (using other methods) [H]. 
We first ported Knuth's Algol program to C++ and quickly found that it was 
prohibitively hard to generate additional numbers given the fact that P3(50) is of 
the order of 10^^. So we decided to parallelize Knuth's program in the following 
way. 

1. Generate all almost topological sequences up to a depth k. 

2. Next, separately run each sequence (to generate the rest of tree) from depth 
(A; + 1) until all sequences of index that contain the initial sequence as 
its first k terms are generated. Here it is important to note that while we 
are counting the numbers of topological sequences, we need to include all 
almost topological sequences since they necessarily appear as sub-sequences 
of topological sequences. 

3. An important observation is that it suffices to run one sequence for every 
given shape since they have identical tree structure after the {k + l)-th 
node. However, it is crucial to note that each topological sequence in a 
given equivalence class does not have the same index. This entails a bit of 
book keeping where one keeps track of the different indices of all topological 
sequences of identical shape. The power of this approach is best illustrated 
by looking at Table [2] where we list the numbers of actual sequences (nodes) 
as well the number of shapes. A naive estimate (based on the reduction of 
the number of runs) shows that run times should go down by an order of 
10^ - 10^ 

This approach has enabled us to extend the Knuth-Mustonen-Rajesh results 
to all integers A^ < 68. The numbers were generated in several steps: A^ = 
52, 55, 62, 68. The results for A^ < 52 we obtained without parallelization. The 
results for A^ < 55 were obtained using parallelization to depth 7 but without 
using equivalence classes and required about 1500 hours of CPU time. The results 
for A^ < 62 were done using parallelization to depth 14 (4167 shapes) and took 



16 



around 30000 hours of CPU time(about a month of runtime). The last set of 
results for < 68 took around 360-ft" hours of runtime (spread over five months). 

We also extended the numbers for four- dimensional partitions of iV < 35 and 
five-dimensional partitions of iV < 30. This was done without any parallelization. 
The complete results are given in appendix O 

5 Conclusion 

We believe that our results show that it is indeed possible to understand the 
asymptotics of higher dimensional partitions. The preliminary nature of our 
results shows that a lot more can and should be done. Our results provide 
a functional form to which results from Monte Carlo simulations, of the kind 
carried out by Mustonen and Rajesh |T4], can be fitted to. However, the errors 
should be better than one part in 10^ or 10^ to be able to fix the sub-leading 
coefficients. We are indeed making preliminary studies to see whether one can 
achieve this. 

Another avenue is to see if there are sub-classes of partitions that can be 
counted i.e., we can provide simple expressions for their generating functions. 
For instance, the analog of conjugation in usual partitions is the permutation 
group, Sd+i, for (i- dimensional partitions. Following Stanley [26], we can organise 
(i-dimensional partitions based on the subgroups of S^+i under which they are 
invariant (see also ^27j). Some of these partitions might have simple generating 
functions. 

One of the proofs of the MacMahon formula for the generating function of 
plane partitions is due to Bender and Knuth |28](see also [29]). It is done by 
considering a bijection between plane partitions and matrices with non-negative 
entries. There is a natural generalization of such matrices into hypermatrices - 
these hypermatrices are counted by MacMahon numbers. It would be interesting 
to contruct a Bender-Knuth type map between solid partitions and hypermatrices 
and study how it fails to be a bijection. This might explain why the asymptotics 
of MacMahon numbers works so well for higher- dimensional partitions. 

Acknowledgments: We would like to thank Arun Chaganty, Prakash Mohan, S. 
Sivaramakrishnan as well the other undergraduate students of IIT Madras' Boltz- 
mann group who provided a lot of inputs to the project on the exact enumeration 
of solid partitions (http://boltzmann.wikidot.com/solid-partitions). We 
thank the High Performance Computing Environment (http : / /hpce . iitm . ac . in) 
at IIT Madras for providing us with a stable platform (the leo and vega super- 
clusters) that made the explicit enumeration of higher- dimensional partitions 
possible. We thank Nicolas Destainville for drawing our attention to ref. |24j . 
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A Asymptotics of the MacMahon numbers 



In this appendix, we work out the asymptotics of the MacMahon numbers using 
a method due to Meinardus ^30j . A nice introduction to this method is found in 
the paper by Lucietti and Rangamani [10]. 

We have seen that the generating function for (i-dimensional MacMahon num- 
bers is given by 

oo oo 

MM = 1 + E ^'^(^) = n - — zj^^ ■ (A-i) 

Inverting this, we obtain: 



n=i „=i(l-g"U.-i 



where g is a complex variable and F is a circle |g| = e < 1 traversed in the 
counterclockwise direction. We shall evaluate the contour integral in ( ]A.2[) by 
writing q = e~* and then taking the limit t 0. This corresponds to the 
contribution to (1A.2P due to the pole at g = 1, which is the dominant contribution. 
The poles of Md{q) occur precisely at all roots of unity, with the sub-dominant 
contributions coming from other roots of unity. 
We have. 



logMd(e- 



f;a„log(l-e-*"), «n=(''^^/)- (A.3) 

n=l ^ ' 



We expand the logarithm inside the sum using its Taylor series and using the 
Mellin representation of i.e.. 



2m 

We obtain 



1 /'7+ioo 

e-^ = — / dsx-'T{s) , 7>0. (A.4) 

7— ioo 



1 /"yi-ioo 

logM,(e-*) = — / dsTis)as + l)D,is)t'\ (A.5) 

where the Dirichlet series Dd{s) defined as 

00 



n=l 



The real constant 7 is chosen to lie to the right of all poles of Dd{s) in the s-plane. 
For d = 3, ttn = n{n + l)/2 and hence the Dirichlet series is 

n=l 
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Hence, D^^s) has simple poles at s = 2, 3 with residue 1/2 at both poles. For 
general d, Dd{s) has poles at s = k, k = 2,3, . . . ,d. Let us denote the residue at 
s = khj Ak. 

Now, we shift the contour in (IA.5I) from Re(s) = 7 to Re(s) = —a, for 
< a < 1. In the process, log Md{q) receives contributions from the poles of the 
integrand that lie between Re(s) = 7 and Re(s) = —a. Hence, we get 

d 

logMrf(e-*) = J] r(A;) ({k + 1)^-'= + D',{0) - 0,(0) logt 



k=2 



+ 7^1 dsT{s)Cis + l)Ddis)t-\ (A.6) 



27ri 

The integral can be shown to go as Hence, we get 

d 

M,(e-*) = exp ( 5^ Ak T{k) ({k + 1) + D'M - D,(0) logt) (l + 0(|t| 

k=2 

Hence, near g = 1, we have 



(A.7) 



1 



(*to-i~i7v 

mJn) = — I dte^-'^^\ (A. 



where (to is taken to close to O'*") 

Gd{t) = 5^ ^ + Ck:=Ak T{k + 1) Cik + 1) 

k=2 

We carry out the integral (]A.8|) using the saddle point method. For this, we 
have to first evaluate t = such that (^^(t*) = 0. That is, 

k=2 ^* 

We next let the integration contour pass through the saddle point for which the 
value of Gdit^) is largest. This happens when is the largest root of (]A.9|) . This 
means t* ^'^^^'^ ~ or equivalently, ~ ^-i/rf+i g^^^ hence, — ?■ as n — )■ 00. 
Hence, the saddle point method indeed gives the value of md{n) for n — )■ 00. 

Now, we solve for t^, from flA.Qp which is a polynomial equation of degree d+1. 
For (i > 3, we do not have a general formula for the roots of the equation. But 
in this case, we indeed have a formula for the largest positive root of flA.Qp . due 
to Lagrange: 

00 

t,{n)= hen-^/^^^+^^ (A. 10) 

e^O mod 
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where 



'' = 1 



d \ d+l 

with (p{y)^[yCky''-'' 



y=0 \ k=l 



Using the above formula, we can compute to to any required order in n and then 
carry out the saddle-point integration (lA.9p . We finally get 



Mn) = pj.^t^) exp {G,iQ + D'M) (l + (C) ) • (A.ll) 

Recall that the dependence on n occurs implicitly, on the right hand side of the 
above equation, through the saddle-point value t^{n). 

A.l Three-dimensional MacMahon numbers 

The asymptotic formula is 

m3{n) ~ const n~^^^^^ exp(G3(n)) , (A. 12) 

wher^ 

Gsin) := ICs^/^n^/^ + Au^'^' ^ {^C,C, - Cl) 
with Ci = 0, C2 = C(3) and C3 = 3^(4). Numerically evaluating, we obtain 

G3(n) ~ 1.78982n=^/^ + 0.333546^^ - 0.0414393n^/^ . (A.13) 

A. 2 Four-dimensional MacMahon numbers 

The asymptotic formula is 

7714(72) ~ const 

2179/3600 exp(G4(r7)) , (A.14) 

where 

r — 1/5^4/5 I C'jn^/^ ^ (5C2C4 - C|) ^2/5 I (^^l ~ 5C2C4C3 + 25CiC|) ^^/^ 

with Ci = 0, C2 = 2C(3)/3, C3 = 3C(4) and C4 = 4C(5). Numerically evaluating, 
we obtain 

Gi{n) ~ 1.66139 n^'^ + 0.460969 n^'^ + 0.0829315 77^/^ - 0.0345152 77^/^ . (A.15) 



^We add a term corresponding to fc = 1 with coefficient Ci in Eq. (jA.9| so that the saddle 
point computation can be carried over for higher-dimensional partitions for which that might 
be the case. 
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A. 3 Five-dimensional MacMahon numbers 

The asymptotic formula is 

m^in) ~ const n-^^^/^^'^ exp(G5(n)) , (A. 16) 

where 

50 ^^2/3 ^2C|/2 

_l_ (2C| — 9C3C5C4 + 27C2CI) ^2/6 

(-QIQ^ + 504C3C5CI - 864^2^1^4 + 432C| (I2C1C5 - Cj)) 



5184C5 



19/6 



with Ci = 0, C2 = K(3)/3, C3 = ^C(4), C4 = 6C(5) and C5 = 5C(6). Numeri- 
cally evaluating, we obtain 

= 1.5737 525874 15873 v^+0.0223817 n^/^-0.0263759 n^/*' . 

(A.17) 



B A rather exact formula for m^{n) 

We will work out the asymptotics of the three-dimensional MacMahon num- 
bers using methods due to Almkvist [2T1[22]. The generating function of three- 
dimensional MacMahon numbers is 

00 00 
M3(x) = JJ(1 - _ ^^3(^) a;- . (B.l) 

n=l n=0 

The integrals are evaluated using the circle method due to Hardy and Ramanujan 
[1]. The coefficients 77^3 (n) are determined from the generating function by the 
formula 

^3(^) = ^ f M,{e'y) e-'^ydy . (B.2) 

Since M^{x) has poles when ever a; is a root of unity, the dominant contributions 
occur in the neighborhood of this point. Setting x = exp{iy), we see that the 
poles occur for all y = 27rh/k with {h, k) = 1 the contribution can be evaluated 
by summing over contributions from such terms. One writes 

00 fc— 1 „ 

m3(n)~^ ]vf3(e'(27rVfc+^)) e-m(2Wfc+^)^y, ^ (B.3) 

k=i h=i -^^h.k 

{h,k)=l 

00 

r-Y^Mn) (B.4) 

k=l 
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where '-)h,k is an arc passing through = 0. We don't give a detailed discussion 
on the choice of the arc but refer the interested reader to [31]. In the second 
hne, we have imphcitly assumed that the integrals and the sum over h have been 
carried out. 

In order to carry out the integral for a particular {h, k), we need to compute 
the Laurent expansion of M^lx) about the point x = exp{2'Kih / k) and then com- 
pute the integral using methods such as the saddle point. For usual partitions, 
this is typically done using modular properties of the Dedekind eta function. 
However, there is no such modular property in this case. The dominant contri- 
bution occurs for /c = 1 (or a: = 1) and we will first consider this contribution. 
Let 

oo oo 

g^ait) := logM3(e-*) = -1 ^ z/(z/ + 1) log(l - e"^*) ^ ^ h^^), (B.5) 

v=l u=l 

where h^dix) := —EiElll log(l — e~^*).The Abel-Plana formula enables us to 
replace the discrete sum over u by the integral: 

93,{t) = h{x) dx-^j^ M^y)^- K-w) ^B_Q^ 

For hr{x) := — x'' log(l — e~^*) , by expanding out the logs and resumming, 
Almkvist has shown that [22] 

, oo 

9r{t) = + C'(-r) - a-r) logt + ki-r -l)]+Yl ^^^"'^^g'^'^^ t- , 

= grit) + gr^'it) , (B.7) 

where in the second line g^^^'it) refers to terms appearing as the sum in the first 
line and gr{t) the remaining terms (within square brackets) up to order t. This 
separation is useful in computing the saddle-point where we will drop the terms 
appearing in gf^'^it) in computing the location of the saddle point. Then, it 
follows that 



gzd{t) = ]^(gi{t)+g2{t)) 



M^{e ) ~ exp 



gi{t)+g2(t) 



(B.8) 



Note that the infinite sum for g2{t) vanishes since C,{—2n) = for n = 1, 2, 3, . . . 
while for gi{t) only terms with even u contribute. In computing the integral in 
Eq. (IR3|), we 

1 r .|[C'{-l)+C'{-2)] roo / ai 2a2 \ 

— / M,{e'^) d^=- / (-,y,)-7el2(-^)^ + 2(-.,,)3 

"'71 1 J-oo 

(B.9) 



22 



where ai = C(3), aa = 2C(4), 7 = C(-l)/2 = -1/24 and ^ = n + Using the 

expansion 



and the integral 



2^1 +^2 z/i ! z/o ! ( -icp) ^''i 



2nJ.J ^ I 5(0 ifa = ^ ^ 



we find that the contribution ignoring the terms in glY^{t) is given by 
0fHn)~exp(l[C'(-l) + C'(-2)]) E 



0^02' ^2.1+3.2-1+7 



2-i+-2zy,!i.2!r(2i/i + 3z/2 + 7) 

(.i,.2)GN*^ 



:=exp(i[C'(-l) + C'(-2)]) L[e,7] , (B.12) 

where we have imphcitly defined the function i^[^,7] in the second fine. In order 
to include the contribution of g^Y^{t), we consider the Taylor expansion (Note 
that Co = 1) 

00 

exp(^/^rW) =E^^-^' ' (B.13) 
i=o 

and carry out the integrations to obtain 



01 (^) = 



00 



n] 



j=0 



:= exp (|[C'(-1) + C'(-2)]) 5^ c, L [n + ^,7 - j] (B.14) 

j=0 

B.l Other poles 

Let us evaluate M3d(e*^) in the neighbourhood of such a point. Put y = 2Txh/k+Lp 
and using a method due to Almkvist(see Theorem 5.1 in |22]), we get 

M^,{e^'-^''-^^) ~ exp (1 [^{-i^r^ + fl(-^V^)-'] + MK'(-I) + k\'{-2)\ 
+ f [s(l, /i, k) + s(2, /i, A;)] - ^C(-l) log(-^M - K(-3)^<^ + ■■■), (B.15) 



23 



where the generahzed Dedekind sums are 



k ik'^t 
s{l, h,k) = -J2 Mj/k) log |2 sinO'/iTT/A;) I + ^ J] B^{j/k) cot{jhn/k) 

fc-i , fc-i 



s{2,h,k) = -Y,B,{3/k){{3h/k)) = --L^cotM(jWfc)cot(j7r/fc) , 

j=l 3=1 

(B.16) 

where Bn{x) are the BernouUi polynomials and 

We illustrate the computation of 0i(n) for n — 60. Below, we quote the result 
after rounding off to the nearest integer and underline the number of correct 
digits. 

0f^(6O) = 11031 748252850258 

<f)f\60) + <f)[^\60) = 1103128 7052778130 
01 (60) = 110312866 33959406 
01 (60) + 02(60) = 11031286641 929870 
m3(60) = 11031286641714044 

Wc observe that (f)f'\60) gets the first five digits right while 0i(6O) makes the 
estimate correct to nine digits while adding 02(60) gets 11 digits right. We need to 
include the contributions of of other zeros i.e., 0fc(n) for A; > 2 to further improve 
the estimate. We anticipate that addition of other terms should eventually lead 
to an exact answer though we have not explicitly verified that it is so. 



C Exact enumeration of higher-dim. partitions 

In this appendix, we provide the results obtained from our exact enumeration 
of three, four and five-dimensional partitions. In all cases, we have gone signifi- 
cantly beyond what is known and we have contributed our results to the Online 
Encyclopedia of Integer Sequences(OEIS) - the precise sequence is listed in the 
table. We believe that it will be significantly harder to add to the numbers of sohd 
partitions as the generation of the last set of numbers took around five months. 
In this case, adding a single number roughly doubles the runtime. There is, how- 
ever, some scope for improvement for the four and five-dimensional partitions as 
the numbers were generated without parallelization. 
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n 




n 




n 







1 


23 


19295226 


46 


8683676638832 


1 


1 


24 


35713454 


47 


14665233966068 


2 


4 


25 


65715094 


48 


24700752691832 


3 


10 


26 


120256653 


49 


41495176877972 


4 


26 


27 


218893580 


50 


69531305679518 


5 


59 


28 


396418699 


51 


116221415325837 


6 


140 


29 


714399381 


52 


193794476658112 


7 


307 


30 


1281403841 


53 


322382365507746 


8 


684 


31 


2287986987 


54 


535056771014674 


9 


1464 


32 


4067428375 


55 


886033384475166 


10 


3122 


33 


7200210523 


56 


1464009339299229 


11 


6500 


34 


12693890803 


57 


2413804282801444 


12 


13426 


35 


22290727268 


58 


3971409682633930 


13 


27248 


36 


38993410516 


59 


6520649543912193 


14 


54804 


37 


67959010130 


60 


10684614225715559 


15 


108802 


38 


118016656268 


61 


17472947006257293 


16 


214071 


39 


204233654229 


62 


28518691093388854 


17 


416849 


40 


352245710866 


63 


46458506464748807 


18 


805124 


41 


605538866862 


64 


75542021868032878 


19 


1541637 


42 


1037668522922 


65 


122606799866017598 


20 


2930329 


43 


1772700955975 


66 


198635761249922839 


21 


5528733 


44 


3019333854177 


67 


321241075686259326 


22 


10362312 


45 


5127694484375 


68 


518619444932991189 



Table 3: Numbers of solid partitions. This is sequence A000293 in the OEIS [32] . 
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n 


P4(n) 


n 


P4(n) 


n 







1 


13 


181975 


25 


2569270050 


1 


1 


14 


425490 


26 


5427963902 


2 


5 


15 


982615 


27 


11404408525 


3 


15 


16 


2245444 


28 


23836421895 


4 


45 


17 


5077090 


29 


49573316740 


5 


120 


18 


11371250 


30 


102610460240 


6 


326 


19 


25235790 


31 


211425606778 


7 


835 


20 


55536870 


32 


433734343316 


8 


2145 


21 


121250185 


33 


886051842960 


9 


5345 


22 


262769080 


34 


1802710594415 


10 


13220 


23 


565502405 


35 


3653256942840 


11 


32068 


24 


1209096875 






12 


76965 


25 


2569270050 







Table 4: Numbers of four-dimensional partitions. This is sequence A000334 
the OEIS [32]. 



n 




n 




n 


P5{n) 





1 


11 


119140 


22 


3923114261 


1 


1 


12 


323946 


23 


9554122089 


2 


6 


13 


869476 


24 


23098084695 


3 


21 


14 


2308071 


25 


55458417125 


4 


71 


15 


6056581 


26 


132293945737 


5 


216 


16 


15724170 


27 


313657570114 


6 


657 


17 


40393693 


28 


739380021561 


7 


1907 


18 


102736274 


29 


1733472734334 


8 


5507 


19 


258790004 


30 


4043288324470 


9 


15522 


20 


645968054 






10 


43352 


21 


1598460229 







Table 5: Numbers of five- dimensional partitions. This is sequence A000390 
the OEIS [32]. 
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